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Abstract 


The impact of Asian dust on cloud radiative forcing during 2003-2006 is studied by using 
the Earth’s Radiant Energy Budget Scanner (CERES) data and the Fu-Liou radiative transfer 
model. Analysis of satellite data shows that the dust aerosol significantly reduced the cloud 
cooling effect at TOA. In dust contaminated cloudy regions, the 4-year mean values of the 
instantaneous shortwave, longwave and net cloud radiative forcing are -138.9, 69.1, and -69.7 
Win' 2 , which are 57.0, 74.2, and 46.3%, respectively, of the corresponding values in more 
pristine cloudy regions. The satellite-retrieved cloud properties are significantly different in the 
dusty regions and can influence the radiative forcing indirectly. The contributions to the cloud 
radiation forcing by the dust direct, indirect and semi-direct effects are estimated using 
combined satellite observations and Fu-Liou model simulation. The 4-year mean value of 
combination of indirect and semi-direct shortwave radiative forcing (SWRF) is 82.2 Wm' 2 , 
which is 78.4% of the total dust effect. The direct effect is only 22.7 Wm' 2 , which is 21.6% of 
the total effect. Because both first and second indirect effects enhance cloud cooling, the 
aerosol-induced cloud wanning is mainly the result of the semi-direct effect of dust. 

Key words: Dust aerosol, cloud, indirect effect, radiation forcing, CERES. 
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1 . Introduction 


It is well known that aerosols have a profound impact on the global and regional 
climate. Determining the radiative effects of aerosols is currently one of the most active 
research areas in climate study (IPCC, 2007) because aerosols can have significant 
impacts on regional and global climate. Aerosols can reflect and absorb the incoming 
solar radiation, which is called the aerosol direct radiative effect. They can increase cloud 
albedo and suppress precipitation by modifying cloud microphysical properties, 
phenomena that are defined as the aerosol first and second indirect radiative effects, 
respectively (Twomey, 1977; Albrecht, 1989). Aerosols can also evaporate low-level 
clouds through so-called aerosol semi-direct effect (Grassl, 1975). Many observational 
and modeling studies have focused on the aerosol ‘direct’ radiative effect, but recent 
research indicates that small changes in cloud microphysical properties (e.g., droplet size, 
phase) could induce large changes in cloud radiative forcing. Therefore, understanding 
and quantifying both the direct and indirect radiative effects are critically important in 
climate research (IPCC, 2007). 

Dust storms that originate in Taklamakan Desert of China and the Gobi Desert of 
Mongolia occur frequently in recent years and become one of important aerosol types in 
East Asia (Huang et al., 2006a). The dust layers associated with these storms often travel 
thousands of kilometers at high altitudes moving from the continent to the open sea near 
Korea and Japan (Haywood et al., 1999; Higurashi and Nakajima, 2002; Takemura et al, 
2002). Moreover, Asian dust particles, which could lead to both cooling by reflecting 
sunlight back to space and warming by absorbing solar and infrared radiation, make the 
radiative impact of Asian aerosols more complicated (Sokolik and Toon, 1999). Thus, it 
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is imperative to examine the aerosol impact on radiative energy budget over East Asia 
where the aerosol sources are presumed to be much different from those in other regions. 

Although there have been numerous studies about the effect of dust aerosols on the 
radiative energy budget in clear skies, few have investigated the radiative effects of Asian 
dust in cloudy conditions. Huang et al. (2006a) studied the impacts of dust storms on 
cloud properties and radiative forcing (RF) using data collected by the Moderate 
Resolution Imaging Spectroradiometer (MODIS) and Clouds and the Earth’s Radiant 
Energy System (CERES) instruments on the Aqua and Terra satellites. They found, on 
average, ice cloud effective particle diameter, optical depth and ice water path of cirrus 
clouds under dust-polluted conditions are 11%, 32.8%, and 42% less, respectively, than 
those derived from ice clouds in dust-free atmospheric environments. Due to the altered 
cloud microphysics, the instantaneous net cloud RF changed from -161.6 W/m 2 for dust- 
free conditions to -118.6 W/m 2 for dust-contaminated clouds. However, the above- 
mentioned results were based on only 10 cases analysis. In this study, we will extend the 
study of Huang et al. (2006a) by considering more cases, and further separate the dust 
direct effect from the aerosol indirect and semi-direct effects on cloud radiative forcing 
using the radiative transfer model along with satellite data analysis. By examining the 
influence of Asian dust aerosols on the radiative forcing both at the TOA and surface, this 
study should lead to a better understanding of the interaction among dust aerosols, clouds 
and radiation. 

2. Data and Methodology 

Four years (June 2003 to June 2006) of CERES Aqua Edition IB Single Scanner 
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Footprint (SSF; see Geier et al., 2001) data are used here. CERES SSF data sets combine 
CERES radiation measurements, MODIS cloud microphysical retrievals (Minnis et al. 
2007) and ancillary meteorology fields to form a comprehensive, high-quality 
compilation of satellite-derived cloud, aerosol, and radiation budget information for 
radiation and climate studies. There are about 140 parameters in the SSF data set. The 
current analysis uses the following SSF parameters including IWP, LWP, cloud top 
effective temperature (T e ), cloud effective droplet radius (Re), ice crystal diameter (De) 
and optical depth, which were derived with the Visible-Infrared-Solar-infrared-Split- 
window Technique (VISST; see Minnis et al., 2004, 2007). Flux measurements from 
CERES are used to estimate the instantaneous RF. The CERES instruments measure 
broadband radiances at the Top-Of-Atmosphere (TOA) in three spectral regions (0. 2-5.0 
pm; 8-14 pm; 5-100 pm) with a spatial resolution of about 20 km at nadir direction. 
These radiances are converted to TOA fluxes with a scene identification algorithm and 
angular distribution models (Loeb et al, 2005). 

To estimate the impact of Asian dust on both cloud radiative forcing and cloud 
microphysical properties, two regions were selected to represent the clouds in different 
environments in the data analysis. CLD (pure cloud) represents the area where clouds 
occurred in a dust-free atmosphere, while COD (cloud- over-dust) denotes the clouds are 
contaminated by dust. The CLD and COD regions are determined based on both MODIS 
images and surface meteorology station observations over China and Mongolia (Huang et 
al., 2006a). The surface stations report dust events in four categories: dust storm, wind- 
blown sand, floating dust, and no-dust. A cloudy region observed by satellite is defined as 
CLD, if the surface observation is no-dust, or as COD, if the surface observation is dust 
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storm, wind-blown sand, floating dust. Table 1 shows information for the 16 elected dust 
cases used in this study. These selected cases occurred during the 2003-2006 time period. 

The concept of cloud radiative forcing is a parameter that has been extensively used 
to study cloud-radiation interaction. Here, we use the outgoing short-wave (SW) and 
long-wave (LW) TOA fluxes to examine the cloud radiative forcing. The cloud radiative 
forcing is defined as the difference in the fluxes between the clear-sky and actual sky 
conditions (Charlock and Ramanathan, 1985; Hartmann et ah, 1986; Ramanathan, 1987; 
Cess and Potter, 1987; Ramanathan et al., 1989.). 

The SW and LW fluxes are defined at the TOA and surface levels as 

rp j-i down 771 up 

1 SW ~ r SW - r SW 

rp t-t down j-y up (i) 

1 LW r LW ^ LW 

and cloud radiative forcing is given as 

✓~r rp rj-i CLR 

SW — 1 SW ~ 1 SW 
rp rpCLR 

LW — 1 LW 1 LW (2) 

C = C +C 

where F^' v " ( F^ v ) and F r d ° wn ( F"^ ) are downward (upward) SW and LW fluxes, 

respectively. The index “CLR” denotes clear-sky scenes without clouds and dust aerosols. 
Thus, the cloud radiative forcing derived in this study for COD conditions includes the 
effects of both clouds and dust aerosols on the radiative energy budget. 


3. Radiation Transfer Model 

The radiation transfer model originally developed by Fu and Liou (1992, 1993) and 
modified by Rose and Charlock (2002) is used to to elucidate the observation results. It is 
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a delta-four stream radiative transfer code with fifteen spectral bands from 0.175 to 4.0 
pm in SW and twelve LW spectral bands between 2850 and 0 cm' 1 . The correlated k- 
distribution method is used to treat the non-gray gaseous absorption due to H 2 0, C0 2 , O 3 , 
N 2 0, and CH 4 (Fu and Liou, 1992). The surface albedo spectral dependencies are taken 
into account using a lookup table based IGBP scene type. The radiation model allows for 
the input of multiple wavelength aerosol optical depth measurements while the aerosol 
wavelength dependent single scatter albedo and asymmetry parameter are determined by 
choosing an assumed aerosol type or external mixture of constituents. 

4. Result Analysis 

Figure 1 shows the histograms of the annual mean instantaneous TOA CRF for dust- 
free cloud (CLD) and dusty cloud (COD) regions for 2003-2006. In each of these 4 years, 
the (absolute) values of SW, LW and Net forcing in COD region are all smaller than 
those derived from CLD areas. The 4-year average values of SW, LW and Net CRF are - 
138.9, 69.1, and -69.7 Wm' 2 , respectively, for COD regions. These values are only 57.0, 
74.2, and 46.3% of those in CLD regions, respectively. Figure 1 indicates that the dust 
reduces the cloud modulation of the radiative energy budget including the cloud albedo 
(cooling) and greenhouse (warming) effects as well as the total cloud cooling effect. 
These results are consistent with our previous study (Huang et ah, 2006a) except that the 
net CRF in the COD regions is 40% lower. It may be due to stronger dust storms during 
2005 and 2006 than during 2003 and 2004. Both previous studies (Huang et ah, 2006a, b) 
suggest that the existence of dust under clouds significantly reduces the cooling by 
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clouds. The reduced cooling due to dust in the COD regions can be considered as a 
warming effect of these aerosols. 

Figures 2 and 3 compare the instantaneous SW, LW and Net CRF at TOA from the 
COD and CLD regions as functions of solar zenith and cloud optical depth, respectively. 
The absolute values of SW and Net CRF in Figs. 2(a,c) and 3(a,c) in COD regions are 
significantly smaller than those in CLD areas. However, for LW CRF, as shown in Figs. 
2(b) and 3(b), the differences between the COD and CLD CRF values are less distinct, 
especially in Fig. 3(b) where the cloud optical depth is smaller than 35. The CRF 
differences between the COD and CLD conditions depend more on the cloud optical 
depth than on the solar zenith angles (Figs. (2) and (3)). 

The above analyses clearly show that cloud radiative forcing is influenced by 
aerosols when dust plumes exist under the cloud. These influences include direct, indirect 
and semi-direct effects of dust aerosols. Figure 4 shows the frequency distributions of ice 
cloud effective diameter (De), water cloud effective radius (Re) and cloud optical depth 
for the COD and CLD regions. The mean values of De and Re for COD regions are 
17.0% and 30.0%, respectively, less than those derived from CLD region. In the dust 
polluted area (COD), the dust may serve as extra condensation nuclei which restrain 
cloud particle growth and lead to the smaller values of De and Re compared to those in 
CLD regions. 

The smaller De/Re associated with the aerosol first and second indirect effects 
(Twomey, 1977; Albrecht, 1989) would enhance the cloud optical depth and cloud solar 
albedo. (Note that the cloud optical depth is proportional to the IWP/LWP but inversely 
proportional to the effective size of cloud particles.) The average COD cloud optical 


8 



depth, 10.1, is, however, 34% smaller than that, 15.6, for CLD areas. The decrease in the 
optical depth may be explained by the evaporation of low clouds due to the dust aerosol 
semi-direct effect that dominates the first and second indirect effects. Such a decrease 
may be also partly attributed to the aerosol radiative effects on surface fluxes and the 
hydrological cycle (Jiang and Feingold, 2006), although their relative contributions to the 
decrease of cloud optical depth remain as an important topic for future investigation. It is 
evident that the dust aerosol can change the cloud micro- and macro-physical properties 
and, thus, the CRF. 

To further elucidate the dust effects, simulations were conducted using the Fu-Liou 
radiative transfer model. In the model calculations, the dust layer was placed below 4 km: 
the aerosol loading decreases exponentially with height and at 4 km it is 1/e of the surface 
value. The type of dust was selected as coarse mode. Vertical profiles of temperature, 
pressure, humidity, and gaseous concentrations correspond to the mid-latitude summer 
(MLS) atmosphere of McClatchey et al. (1971). The cloud was placed between 3.5 and 
5.5 km over a desert surface. 

The direct effect of dust aerosols consists of scattering and absorption of incoming 
solar radiation. Because dust is an absorbing aerosol, it can exert different influences on 
the SW and LW fluxes at both the TOA and the surface. As dust optical depth (x) 
increases (Fig. 5), the cooling related to the SW and Net CRF decreases at the TOA, 
while it increases significantly at the surface. These differences imply that the dust 
aerosols decrease the cloud cooling effects at TOA and increase the cloud cooling effects 
at surface. The reduced/enhanced cloud cooling effects are actually the direct 
warming/cooling effects due to dust aerosols. From Fig. 5, it is clear that when the 
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aerosol optical depth increases by 0.1, the value of instantaneous SW CRF will change 
from 2.7 to 15.7 Wm' 2 at TOA and -1.7 to -29.7 Wnf 2 at surface, depending on the solar 
zenith angles. When dust aerosols occur under clouds, their SW and Net warming effects 
at TOA are due to the dust absorption, while both the absorption and scattering by dust 
prevent solar energy from reaching the surface and cause a cooling there. The difference 
in LW CRF between dust and dust-free cloudy skies is always positive because of the 
greenhouse effect of dust aerosols. But it is small because the long wave is relatively 
insensitive to the presence of dust. The impact of these aerosols on LW radiation at the 
TOA is smaller than that at the surface because the clouds above the aerosol layer mask 
any dust effect. 

Figure 6 shows the instantaneous SW, LW and Net cloud radiative forcings at the 
TOA and surface as functions of cloud optical depth over desert surface. When the cloud 
optical depth changes, the dust still has a warming effect at the TOA and cools the 
surface. These simulation results suggest that as an absorbing aerosol, dust has the direct 
effect of diminishing cloud cooling at the TOA, so, i.e., dust has a warming direct effect 
at the TOA. 

Satellite observations suggest that the total dust aerosol radiative forcing, which 
includes direct and, through cloud property changes, semi-direct and indirect effects (i.e., 
the difference in CRF between COD and CLD regions), is positive at TOA. However, it 
is difficult, using only satellite observations, to separate the dust aerosol direct effect 
from those caused by altered cloud properties. To estimate the separate effects, we 
propose a simple method that uses both satellite observations and the model calculations. 
Figure 7 compares the CERES-measured instantaneous SW CRF values with the Fu-Liou 
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model simulation results for the CLD regions. The comparison shows good agreement 
between CERES observations and the Fu-Liou model simulation, indicating that the Fu- 
Liou model can reproduce the TOA radiation flux for dust-free cloud regions reasonably 
well. This result suggests that the contribution of dust aerosol direct and indirect/semi- 
direct effects can be estimated by using the differences between the satellite observations 
and the model simulations. The procedure for estimating the contribution of dust indirect 
and semi-direct effect is as follows: 

1. Retrieve the dusty and dust-free cloud properties (i.e. cloud optical depth, Re, De, 
cloud top pressure) through satellite measurements. 

2. Simulate the dusty and dust-free cloud RF using satellite-retrieved cloud 
properties as input parameters. (Note that the dust optical depth was set to be 
zero in both simulations, although the cloud properties were changed by dust 
aerosols in the COD regions). 

3. The total dust aerosol radiative forcing including direct and semi-direct and 

indirect effects is 

ARF to tai = Fceres(COD) - Fceres(CLD) 

where Fceres(COD) and Fceres(CLD) are the radiative energy budget at the 
TOA from CERES observations under COD and CLD conditions, respectively. 

4. The combination of dust aerosol indirect and semi-direct forcing related to 
cloud property changes is 

ARFindirect+semi — F mo del(COD) — F m0( Jel(CLD) 
where F mo d e i is the radiative energy budget at the TOA from the Fu-Liou model 
computed using satellite-retrieved cloud properties, as in steps 1 and 2. 
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5. The direct radiative forcing due to dust aerosols is 

ARF direct — ARF total “ ARF indirect+semi 

The averaged dust aerosol direct and combination of indirect and semi-direct effect on 
instantaneous SW RF at TOA for the selected cases are shown in Table 2. The averaged 
direct, and combined indirect and semi-direct instantaneous SW RF are 22.7 and 82.2 
Wm' 2 , respectively; which correspond to 21.6 and 78.4% of the total RF value. The dust 
aerosol direct effect of 22.7 Wm' 2 as derived from the satellite-measured radiation and 
cloud properties here is within the range of those from the model as shown in Figs. 5 and 
6 (where the dust aerosol direct radiative forcing is the difference between the 
dashed/dotted line and solid line). The percentages of direct and combined indirect and 
semi-direct RF at TOA for the 16 selected cases are shown in Figure 8. It is obvious from 
Table 2 and Fig. 8 that the combination of dust indirect and semi-direct effects is the 
primary contributor to the dusty cloud warming effect. 

5. Conclusion and Discussion 

Asian dust aerosols have important influences on the radiative energy budget 
through both the direct effect and effects on the cloud microphysical properties. Most 
previous analyses of aerosol radiative forcing were limited to clear sky only (direct 
forcing). Recent studies have found that dust plays an important role in modulating 
climate through the indirect aerosol effect on cloud properties (Sassen et al, 2003). Our 
analyses of satellite observations show that the absolute values of instantaneous TOA 
CRF in COD regions are significantly smaller than those derived from CLD regions. This 
implies that when dust exists under a cloud it has a warming influence at TOA. It is 
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clearly evident from analysis of satellite data that dust aerosols can change cloud 
microphysical properties, such as water droplet radius or effective ice particle diameter 
and cloud optical depth. These properties consistently yield smaller values when derived 
from COD regions as compared to those derived over CLD regions. The results presented 
here are consistent with our previous study (Huang et al., 2006a). The warming effects of 
dust aerosols at TOA through dust aerosol direct effects in COD regions were also 
indicated by Fu-Liou model simulation. 

Although the impact of indirect/semi-direct effect of dust aerosols can be seen from 
the satellite-retrieved cloud microphysical properties, it is difficult to separate those 
effects on the radiative energy budget from their direct effect using satellite observations 
alone. This study has proposed a new approach by combining observation and simulation 
results to quantify the indirect and semi-direct effects of dust aerosols as well as their 
direct effect. Satellite observations along with the radiative transfer calculations show 
that the dust aerosol indirect and semi-direct effects dominate the dust aerosol total RF. 
The four-year (2002-2006) combined average indirect and semi-direct instantaneous 
shortwave radiative forcing (SWRF) is 82.2 Win' 2 , a value that accounts for 78.4% of the 
total RF value. The dust aerosol direct effect of 22.7 Wm’ 2 as derived from satellite 
radiation and cloud observations here is within the range of those determined from the 
model simulations. Because both the first and second indirect effect cause cooling at 
TOA, the warming effect of dust aerosols in cloud-over-dust systems is due to the semi- 
direct effect of dust aerosols. 

The term semi-direct effect was introduced to describe the impact of absorbing 
aerosols on clouds (e.g., Hansen et al. 1997). A series of experiments with a simple GCM 
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by Hansen et al. (1997) showed that increased shortwave absorption could reduce relative 
humidity and, then, decrease cloud cover. Similar results were later obtained by Cook and 
Highwood (2003) with a more sophisticated GCM. When the aerosol-induced cloud 
feedbacks were included, the GCM predicted wanning increased from 2.5K to 2.9K. The 
additional warming was related to decreases in the fractional coverage of large-scale 
clouds, particularly at mid and high latitudes. The semi-direct aerosol effect was also 
investigated by Ackerman et al. (2000) using a Large-Eddy-Simulations (LES) based on 
observations from the Indian Ocean Experiment (INDOEX) during 1998-99. They found 
that absorbing aerosols reduced the relative humidity in the boundary layer and caused a 
5 - 10 % reduction in cumulus cloud fraction. Lohmann and Feichter (2001) provided the 
first assessment of the global annual mean aerosol semi-direct forcing using the European 
Center HAMburg 4 General Circulation Model (ECHAM4 GCM). The strong warming 
influence of absorbing aerosols was also emphasized by Jacobson (2002), who found a 
decrease of global averaged column liquid water and ice associated with the semi-direct 
effect. Recently, Huang et al. (2006b) studied the semi-direct effects of dust aerosols over 
eastern Asia by using 2 years of CERES and ISCCP data. They found that the semi-direct 
effect might play an important role in cloud development over arid and semi-arid areas of 
East Asia and contribute to the reduction of precipitation. 

Although this paper is only a preliminary study of Asian dust effects in cloudy 
conditions, it sets the stage for more in-depth analyses by providing a method to separate 
the various types of aerosol effects. It has also demonstrated how important these non- 
direct effects are in areas with heavy dust loadings. The use of other datasets in 
combination with the model and CERES data will further enhance the quantification of 
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the different effects. For example, the Cloud- Aerosol Lidar with Orthogonal Polarization 
(CALIOPE) on the CALIPSO satellite, launched in April 2006, offers an unprecedented 
opportunity to improve the understanding of aerosol forcing by giving accurate and 
detailed aerosol vertical profdes. Through the use of these data it should be possible to 
quantify the vertical structure of dust aerosol radiative forcing. Future analyses will bring 
together both the active and passive datasets with the model calculations to understand 
the effects and provide datasets that will be valuable for validating and improving the 
GCM estimates of aerosol radiative forcings. 
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Table 1 . Sixteen Dust Case Images Used in This Study 


Case Image 

Date 

Lat. ° N 

Lon. ° E 

1 

2003/03/26 

44.0-47.5 

123.5-124.5 

2 

2003/04/09 

36.5-40.5 

78.5-84.0 

3 

2004/03/09 

44.0-46.0 

104.0-113.0 

4. 

2004/03/27 

43.0-47.0 

111.5-115.5 

5 

2005/03/09 

39.5-45.5 

108.5-111.5 

6 

2005/04/06 

48.5-50.5 

112.5-118.5 

7 

2005/04/08 

40.0-41.5 

102.5-105.5 

8 

2005/04/09 

43.5-46.5 

116.5-118.5 

9 

2005/04/13 

45.5-49.5 

112.0-117.0 

10 

2005/04/20 

36.5-37.5 

116.5-119.0 

11 

2006/02/17 

43.0-48.0 

113.0-115.0 

12 

2006/02/21 

41.5-44.5 

123.0-125.5 

13 

2006/04/17 

36.5-40.5 

100.0-106.5 

14 

2006/04/18 

41.0-48.0 

116.0-122.0 

15 

2006/04/22 

42.0-48.0 

116.0-117.0 

16 

2006/04/28 

43.0-45.0 

104.5-111.5 
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Table 2. Comparison of the Averaged instantaneous RF effect at TOA 




Mean SWRF 
Wm' 2 

Standard 

Deviation 

Mean Percentage 
(direct / total) 

Standard 

Deviation 

Direct 

Indirect 

effect 

effect 

22.7 

82.2 

44.45 

41.15 

19.21% 

9.80% 
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Fig. 1. Annual mean instantaneous TOA dusty (COD) and dust-free (CLD) cloud 
radiative forcing for (a) SW forcing, (b) LW forcing and (c) Net forcing. 
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Fig. 2. Comparison of TOA dusty (COD) and dust-free (CLD) instantaneous cloud 
radiative forcing as a function of solar zenith angle for (a) SW forcing, (b) LW forcing , 
and (c) Net forcing . 
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Fig. 3. Same as Fig. 2 but as a function of cloud optical depth. 
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Fig. 4. Comparison of dusty and dust-free cloud properties for (a) ice cloud De, (b) water 
cloud Re and (c) cloud optical depth. The histogram intervals are 5 pm for Re, 20 pm for 
De and 10 for optical depth. 
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Fig. 5. Model simulated dusty instantaneous cloud radiative forcing as a function of solar 
zenith angle for (a) TOA SW forcing, (b) surface SW forcing, (c) TOA LW forcing, (d) 
surface LW forcing, (e) TOA net forcing, and (f) surface net forcing. The amount of 
aerosol loading is indicated on the plots, and cloud optical depth is 16. The aerosol type is 
Coarse Mode. The surface type is desert. 
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Fig. 6. Same as Fig. 5 but as a function of cloud optical depth. 


28 



Fig. 7. Comparison of the instantaneous dust-free cloud CRF value (in CLD region) of 
CERES measurements with Fu-Liou model simulation. 
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Cnse Number 

Fig. 8. The mean values of direct and combination of indirect and semi-direct 
instantaneous SWRF at TOA for 16 cases. The numbers represent the percentage of 
direct (direct/total). 
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